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Measurement of the Inclusive Charm Cross Section at 4.03 GeV and 4.14 GeV 
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The cross section for charmed meson production at y/s — 4.03 and 4.14 GeV has been measured with the 
Beijing Spectrometer. The measurement was made using 22.3 p6~^ of e^e~ data collected at 4.03 GeV and 
1.5 ph~^ of e+e~ data collected at 4.14 GeV. Inclusive observed cross sections for the production of charged 
and neutral D mesons and momentum spectra are presented. Observed cross sections were radiatively 
corrected to obtain tree level cross sections. Measurements of the total hadronic cross section are obtained 
from the charmed meson cross section and an extrapolation of results from below the charm threshold. 

I. INTRODUCTION is the charm threshold region where the hadronic 

cross section has only been measured with an accu- 
The hadronic cross section of e+e" at aU energies ^acy of 15 - 20% [0. Traditionally, Ghadron is mea- 
ls needed to calculate the effects of vacuum polariza- sured by counting hadronic events. This method re- 
tion on parameters of the Standard Model. The en- quires a detailed understanding of trigger conditions, 
ergy region which contributes the largest uncertainty t^e efficiency of hadronic event selection criteria, and 



a subtraction of two photon events and other back- 
grounds. An alternative is measuring (Tcharm and 
adding this to an extrapolation of the 'Ju,d,s contri- 
bution from the region below charm threshold. The 
charmed mesons used in this study are the Z?° and 
_D+. The Ds cross sections are taken from earlier 
works. There is no evidence for continuum charmo- 
nium production [0. The charm counting method 
is intrinsically less sensitive to trigger conditions, 
beam-related backgrounds, and two photon back- 
grounds due to the distinctive topology of charmed 
meson events. 



II. DATA SELECTION 

The data used for this analysis were accumulated 
with the Beijing Spectrometer [^; the total inte- 
grated luminosity at 4.03 (4.14) GeV was 22.3 (1.5) 
pb~^. Candidate tracks were required to have a good 
track fit passing within 1.5 cm of the collision point 
in R and 15 cm in z, and satisfying | cos6'| < 0.85. 
Particle identification was provided by an array of 
time of flight scintillation counters (TOF) and spe- 
cific ionization measurements in the drift chamber 
used for charged particle tracking (dE/dx). For pi- 
ous, consistency {CL{n) > 0.1%) and loose electron 
rejection {Lt,/{L.,^+L(,) > 0.2) were required, where 
Lx is the likelihood for hypothesis X. Kaon identi- 
fication required consistency {CL{K) > 0.1%), pion 
rejection {Lk/{Lk + ^tt) > 0.5) and loose electron 
rejection {Lj^/{Lj^ + L^) > 0.2). Multiple counting 
of D" candidates was removed by positively identify- 
ing pions using the selection {Lk/{Lk + L.^) < 0.5). 
Muons are rejected using a momentum dependent 
criteria based on track penetration into the muon 
detector. 



III. D° AND D+ SIGNAL 

Inclusive K^t:'^ and K~Tr~^T:^ invariant mass dis- 
tributions are shown in Figs. H and for 4.03 and 
4.14 GeV, respectively. (Here and throughout this 
paper, reference to a state also implies its charge con- 
jugate state.) Each histogram was fit to a function 
which included a Gaussian signal plus a background 
function. For the Ktt distributions (Figs. ^ and||a), 
the background function consisted of a Gaussian 
centered at 1.60 GeV to account for the contribu- 
tion to the Ktt spectrum from D" -^ Ktttt^ decays 
plus a third order polynomial background. For the 
Kmr distributions (Figs. Illb and 0b), the background 



function consisted of a third order polynomial. The 
number of signal events in each mode is given in 
Table |. 

The momentum spectra of D candidates with a 
mass between 1.81 and 1.91 GeV is presented in 
Fig. H for ^/s = 4.03 GeV. The three momentum 
regions near 0.15, 0.55, and 0.75 GeV correspond 
to D*D , D*D, and DD production, respectively. 
The spectrum for y/s = 4.14 GeV is shown in Fig. 0. 
The momentum regions near 0.45 and 0.70 and 0.90 
GeV correspond to D*D , D*D and DD produc- 
tion, respectively. The shapes of the D*D spec- 
trum and part of the D* D spectrum are broadened 
due to Doppler-smeared D mesons coming from D* 
decays. In addition, a small low momentum tail on 
each structure is expected due to initial state radi- 
ation. The background shape under the momentum 
spectrum is not flat, making a direct subtraction dif- 
ficult. 



IV. CROSS SECTION 

If the reconstruction efficiency for D mesons were 
constant with respect to momentum, the observed 
cross section could be determined using 



DX) = 



A^(signal) 
eBC ' 



(1) 



where Af(signal) is the number of signal events, e 
is the efficiency, B is the branching fraction of the 
D meson to a decay mode and C is the luminos- 
ity. However, Monte Carlo studies show some mo- 
mentum dependence to the reconstruction efficiency 
(Fig. |5|). In order to measure the cross section, the 
momentum spectrum for D^ and D^ mesons from 
50 to 850 MeV was divided into 20 (40) MeV slices 
for 4.03 (4.14) GeV data. For each momentum slice, 
the invariant mass distribution was fit with a Gaus- 
sian plus a polynomial background. The central 
value of the Gaussian was fixed at the nominal D 
mass; the width was fixed to a momentum depen- 
dent value determined by a fit to a coarser slicing 
of the data (Fig. ^. The differential cross section 
with respect to momentum for this data is shown in 
Figs. and | for 4.03 and 4.14 GeV, respectively. 
The cross section times branching fractions of D° 
and D^ mesons and cross sections calculated using 
the branching fractions of ref. |4| are shown in Ta- 
ble O. The a ■ B values from this measurement are 
compatible with previous measurements by Mark I 
II 1) and Mark II § as shown in Fig. |. 



CORRECTIONS FOR INITIAL STATE 
RADIATION 



gDo{Snom) 






(6) 



The tree level cross sections for charm at y/s ~ 
4.03 GeV and y/s = AAA GeV were obtained by cor- 
recting the observed cross section for the effects of 
initial state radiation (ISR). The ISR correction is 
dependent on the cross section for all energies less 
than the nominal energy. Since these measurements 
were performed only at two energies, some theoret- 
ical modeling of the cross section distribution was 
required. Two different theoretical predictions for 
these cross sections were used in this analysis, the 
Coupled Channel Model l^, and a P-wave phase 
space formalism. The models provide predictions 
for cr_B,i(seff), the tree level (Born) cross section as 
a function of the effective center of mass energy 

squared for production mode i, where i = D'^D , 
D+D-, D*^Tf, D*+D-, D*°D*°, D*+D*-. The 
effective center of mass energy squared is given by 



Scff = Sn 



.(1-fc), 



(2) 



where k ■ E'bcam is the energy of radiated photons and 
Snom is the nominal center of mass energy squared. 
The tree level cross sections are convoluted with a 
sampling function f{k, Snom) that represents a first 
order calculation of the effective luminosity for ra- 
diated photons in a two body radiation model [BJ, 
giving the observed cross section at Snom: 



/o dk ■ f{k, 



i)o-B,j(Soff)(l + (5vp(scff))- (3) 



The vacuum polarization correction (1 + Syp) in- 
cludes both leptonic and hadronic terms. It varies 
from charm threshold to 4.14 GeV by less than ±2%. 
It is treated as a constant with the value 



(l + ,5vp) ==1.047 ±0.024 



(4) 



and moved outside the integrand. With this sim- 
plification, the ISR correction, gi{snom)i is defined 

by 



9i\^nom) — 



•^obSj-i v-^noinj 



f^_B,i('5noin)(l + ^Vp) 



(5) 



The D*° and Z)*+ branching fractions are used to 
calculate Njjoi and Nu+ j, the mean number of D'^ 
and D~^ per production mode i event. The Njjoi and 
Njj+i values are used to weight the ISR correction 
for each mode to obtain the ISR correction averaged 
over all production modes: 



gD+[ST. 



J2i9iiSnom)ND+,^ 



E^N 



D+, 



(7) 



This procedure results in the corrections shown 
in Fig. |l^. Since neither method models the data 
precisely, and the two models vary differently with 
energy, a systematic uncertainty is assigned to be 
one half the rms difference between the two models 
over the energy range 3.9 GeV to 4.2 GeV excluding 
the region from 4.021 to 4.027 where the Coupled 
Channel Model DD cross section is tiny, causing the 
ISR correction to diverge. The corrections for initial 
state radiation at v^ ~ 4.03 GeV and y/s = 4.14 
GeV are: 



4.03 GeV: goo =0.67 ±0.05 
gD+ =0.73 ±0.05 

4.14 GeV: 51^0 =0.83 ±0.06 
gD+ =0.84 ±0.06. 



(8) 

(9) 

(10) 

(11) 



The ISR correction for Dg mesons was calculated 
using the Coupled Channel Model and the same P- 
wave phase space formalism. Figure O shows a pre- 
diction for the tree level cross section of Ds and D* , 
the observed cross section of Ds and D*, the frac- 
tion of Dg mesons from direct production, and the 
fraction of Dg mesons from Z?* decays. Figure [lO| c 
shows the ISR correction for Dg mesons. The ISR 
contribution to the systematic error for Dg produc- 
tion is taken as one half the rms difference between 
the two models: 



4.03 GeV: go^. 


= 0.73 ±0.04 


(12) 


4.14 GeV: .g^. 


= 0.78 ±0.05. 


(13) 



The ISR and vacuum polarization corrections are 
applied to the observed D'~' and D^ cross sections 
found in Table O to obtain the tree level cross section 
for ZJO and D+ at 4.03 GeV and 4.14 GeV as shown 
below: 



4.03 GeV: aoo + cr-^o 



4.14 GeV: cr/jn 

0'D+ - 



(^D- 



= 19.9 ±0.6 ±2.3 nb (14) 
= 6.5±0.2±0.8nb (15) 



^fj-po =9.3±2.1±l.lnb (16) 
ao- = 1.9±0.9±0.2nb. (17) 



Systematic uncertainties are treated in Section VI. 

The BES observed Dg cross section at 4.03 GeV 
is crn+ + an- = 0.62 ± 0.12 ± 0.20 nb fo|. After 



'D- 



applying the ISR and vacuum polarization correc- of these systematic uncertainties are shown in Ta- 
tions, the tree level Z?, cross section is a 



D7 



ble IV 



0.81 ± 0.16 ± 0.27 nb. The Mark III observations of 
Ds at 4.14 GeV give cr^+ +ct^- = 1.34±0.32±0.34 
nb ||ll[. After correction, the tree level value is 
cr^+ + cr^- = 1.64 ± 0.39 ± 0.42 nb. 



VI. SYSTEMATICS 

Several systematic checks were performed. The 
numbers of signal events in the distributions shown 
in Figs. 0, and were compared to the sum of signal 
events from each momentum slice as shown in Ta- 
ble p|. Good agreement for the 4.03 GeV data val- 
idates the slicing technique. For the 4.14 GeV data 
set, which is much smaller, the agreement is poorer. 
This could be due to statistical fluctuations. 

The analysis was repeated using wrong sign com- 
binations (iir+7r+, i^+7r+7r~) to explore systematic 
bias from the slicing and fitting procedure. As ex- 
pected, the structures that are so evident in the right 
sign spectra are absent in the wrong sign spectra. 
(Fig. n2). There is, however, a small excess when 
the K^tt^t:^ spectrum is integrated: 

(o-^o + <7-o)ws = 0.29 ± 0.21 nb (18) 

i(7D++(^D-)ws = 0-7±0.2iih. (19) 

Both the data and the charged MC show a small ex- 
cess that was not present in neutral decays. A much 
larger MC study would be required to prove whether 
the source is procedural, or if there is a small feed 
down from right sign channels into the wrong sign 
analysis. If the source of the excess is procedural, 
the MC efficiency calculation should correct for the 
effect in the data. In either case, no correction is 
required. 

Systematic errors arising from the choice of pa- 
rameters were evaluated by repeating the analyses 
using different bin sizes, fitting ranges, and back- 
ground function shapes. Electron particle identifi- 
cation dominates the integrated luminosity uncer- 
tainty as determined from wide angle Bhabha scat- 
tering events. The uncertainty is evaluated by com- 
paring samples selected independently using dE/dx 
and the barrel shower counter. In addition there 
are systematic errors due to the uncertainties in the 
Monte Carlo-determined reconstruction efficiency, 
errors in the charmed meson branching fractions, 
and the uncertainties in the evaluation of the ISR 
correction scheme as discussed above. Magnitudes 



Sources of systematic uncertainty are segregated 
into components that are common or independent 
for Z?", Z?+, and Ds measurements. The common 
components are the integrated luminosity measure- 
ment, the ISR correction, the vacuum polarization 
correction, and a portion of the D branching fraction 
uncertainties. Since the absolute branching fraction 
scale for _D+ mesons depends on the D^ branching 
fraction scale, the total percentage uncertainty for 
D^ branching fraction (6.7%) is split into a common 
component that matches the percentage uncertainty 
for the D" branching fraction (2.3%) and an inde- 
pendent component (6.2%). All other systematic 
uncertainties are treated as independent and added 
in quadrature. Values are found in Table p^ . 

The total observed D^ and £>+ cross sections are 
shown in Table |l|. Tree level 13° and Z3+ cross sec- 
tions are shown in Table [V[ 



VII. TOTAL INCLUSIVE CHARM CROSS 
SECTION 

Since all D mesons are produced in pairs, the tree 
level non-strange D cross sections are: 

4.03 GeV: a{DDX)= 13.2 ± 0.3 ± 1.4 nb (20) 
4.14 GeV: cr(L>UA)= 5.6 ± 1.1 ± 0.6nb. (21) 

Adding the tree level D cross sections to the tree 
level Ds cross sections gives the total tree level 
charm cross section: 



4.03 GeV: cichar 
4.14 GeV: fJch; 



13.6 ± 0.3 ± 1.5 nb (22) 

arm- 6.4 ± 1.2 ± 0.7 nb. (23) 



These results are compared to Coupled Channel 
model predictions in Table 0. 



VIII. MEASUREMENT OF Rd AND R 

A measurement of Rb is obtained by dividing 2 x 
o'charm by the QED prediction for the tree level muon 
pair cross section 

^ _, 86.8 nb 



a{e e 



+ c- 



Ai M J = 



s (GeV)2 



(24) 



givmg: 

4.03 GeV: i^i, = 5.10 ± 0.12 ± 0.55 (25) 



4.14 GeV: i?i) = 2.53 ± 0.46 ± 0.27. 



(26) 



The contribution to R in the charm threshold region 
from the hght quarks, Ruds is estimated to be 2.5 ± 
0.25 JTJ. This value was compiled from an average 
of measurements of R below charm threshold |1J] 
p3[ p^ . The theoretical expectation is that Ruds is 
approximately independent of center of mass energy 
in this region. The value of R is evaluated using 
R = Rd/2 + Ruds giving: 

4.03 GeV: i? = 5.05 ± 0.06 ± 0.37 (27) 

4.14 GeV: i? = 3.76 ± 0.23 ± 0.28. (28) 

This measurement is more precise, but compatible 
with previous R measurements |T3] |lj] using the 
total cross section method shown in Figures [l^ a and 
Ob and a previous measurement M employing a 
similar Rjj technique shown in Figure [l3| c. Charm- 
counting complements direct R measurements since 
the two methods feature different systematic uncer- 
tainties. 
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Mode 



Sc.m.(GeV) 



Events 



D+ -^ R-TT+TV^ 

D" ^ K-7V+ 

D+ -^ K-TV+TV^ 



4.03 
4.03 
4.14 
4.14 



4174 ± 163 

2341 ± 138 

249 ± 36 

126 ±31 



TABLE I. Number of inclusive signal events for each mode. 



4.03 GeV 



^ = 4.14GeV 



{ao+ +aD-)- B(D+ -^ K--K+-K+) 



0.537 ± 0.015 ± 0.047 nb 

0.449 ± 0.017 ± 0.036 nb 

13.9 ± 0.4 ± 1.3 nb 

5.0 ±0.2 ±0.5 nb 



0.31 ±0.07 ±0.03 nb 

0.15 ±0.07 ±0.01 nb 

8.1±1.8±0.7nb 

1.7 ±0.8 ±0.2 nb 



TABLE II. The observed cross section times branching fraction of Z)'^ and D^ and the observed cross 
section of D" and D+ at 4.03 GeV and 4.14 GeV. 





Inclusive D 


= 4.03 GeV 

Sum of Shces 


Inclusive D 


-- 4.14 GeV 

Sum of Slices 


D'' 


4174 ±163 


4232 ±114 


249 ± 36 


210 ± 46 


D+ 


2341 ± 138 


2395 ± 94 


126 ±31 


66 ±33 



TABLE III. Comparison of the number of events found from the inclusive mass fit and the sum of the 
events from invariant mass fits for the momentum slices. 



Common 


4.03 D" 


4.03 D+ 


4.14 D° 


4.14 D+ 


4.03 Ds 


4.14 Ds 


Luminosity 


5% 


5% 


5% 


5% 


5% 




D Branching Fraction 


2.3% 


2.3% 


2.3% 


2.3% 






ISR 


7% 


7% 


7% 


7% 


6% 


6% 


Vacuum Polarization 


2.2% 


2.2% 


2.2% 


2.2% 


2.2% 


2.2% 


Independent 


4.03 D° 


4.03 D+ 


4.14 i3° 


4.14 D+ 


4.03 Ds 


4.14 Ds 


D Branching Fraction 




6.2% 




6.2% 






MC statistics 


5% 


4% 


5% 


4% 






Fit parameters 


5% 


5% 


5% 


5% 






Previous measurements 










31.9% 


15.0% 



TABLE IV. Common and independent sources of systematic uncertainty. When summing cross sections 
for difi'erent D species, common uncertainties are added linearly and independent uncertainties are added 
quadratically. 



0"do 


+ 


^Z^ 


o"d+ 


+ 


O-D- 


^Dt 


+ 


""d- 



\/s = 4.03 GeV Experiment Coupled Channel Model 

19.9 ± 0.6 ± 2.3 nb 18.2 nb 

6.5±0.2±0.8 nb 6.0 nb 

0.81 ± 0.16 ± 0.27 nb 1.61 nb 

acharm 13.6 ± 0.3 ± 1.5 nb 12.9 nb 

y/s — 4.14 GeV Experiment Coupled Channel Model 

aoo+f^o" 9.3±2.1±1.1 nb 15.1 nb 

oo++(Td- 1.9±0.9±0.2 nb 4.5 nb 

'^D++'^D- 1.64 ± 0.39 ± 0.42 nb 1.85 nb 

(Tcharm 6.4 ± 1.2 ± 0.7 nb 10.7 nb 

TABLE V. Comparison of tree level cross section measurements with predictions of the Coupled Chan- 
nel Model. 
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FIG. 1. The invariant mass of a) Kn and b) 
Kmv tags with a momentum between and 1 GeV 
at Ec.m. = 4.03 GeV. 
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FIG. 3. The momentum of a) Ktv and b) Kmr 
tags with a mass between 1.81 and 1.91 GeV at 
Ec.m. = 4.03 GeV. 
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FIG. 2. The invariant mass of a) Ktv and b) 
Knn tags with a momentum between and 1 GeV 
at Ec.m. = 4.14 GeV. 
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FIG. 4. The momentum of a) Ktv and b) Knn 
tags with a mass between 1.82 and 1.92 GeV at 
Ec.m. = 4.14 GeV. 




FIG. 5. The momentum dependence of the re- 
construction efficiency. 
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FIG. 7. The differential production cross sec- 
tion of a) e+e- -^ D°X and b) e+e" ^ D+X 
versus the momentum of the reconstructed D at 
Ec.m. = 4.03 GeV. 
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FIG. 6. The momentum dependence of the mass 
resolution. 
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FIG. 8. The differential production cross sec- 
tion of a) e+e^ -^ D'^X and b) e+e" -^ D+X 
versus the momentum of the reconstructed D at 
E,.„,. = 4.14 GeV. 
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FIG. 11. a) The tree level cross section for Ds 
mesons, b) The tree level cross section for Dl 

c) The observed cross section for Ds 

d) The observed cross section for D* 
The solid curve was obtained using the 

Coupled Channel model. The dashed line was ob- 
tained using the p-wave phase space formalism. 
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